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1. Introduction 

Calculating the spectrum of onium states is a significant challenge for lattice gauge theory. 
A number of levels can be studied for both charm and bottom quarks. However, dealing with 
heavy quarks requires special care ^ ||]. Using improved staggered sea quarks p|], it is possible 
to reproduce many of the most important features of the spectrum [Q], which had not been done in 
the quenched approximation. This paper updates our work presented at Lattice 2003 



2. Calculational Details 

Ensembles for three lattice spacings were provided by the MILC Collaboration ||]]: a^O. 18 fm 
("extra-coarse"), a ~ 0.12 fm ("coarse"), and ass 0.086 fm ("fine"). (See Table [I].) For the extra 
coarse am q I am s = 0.6, 0.4, 0.2 and 0.1; for the coarse lattice, we also have 0.14, but we have 
only analyzed two values 0.4 and 0.2 for the fine lattice. From 400 to 600 configurations have 
been analyzed in most ensembles. The most notable exception is the coarse ensemble with am q I 
am s = 0.1. For each of the lattice spacings, the scale of each ensemble with different sea quark 
masses was kept approximately fixed using the length r\ [Q, ||] from the static quark potential. The 
absolute scale from the T 2S-1S splitting was determined by the HPQCD/UKQCD group ^ on 
most of our ensembles implying r\ = 0.318(7) fm. 



am q 1 am s 


10/g 2 


size 


volume 


config. 


a (fm) 


0.0492 / 0.082 


6.503 


16 3 x48 


(2.8 fm) 3 


401 


0.178 


0.0328 / 0.082 


6.485 


16 3 x48 


(2.8 fm) 3 


331 


0.177 


0.0164 / 0.082 


6.467 


16 3 x48 


(2.8 fm) 3 


645 


0.176 


0.0082 / 0.082 


6.458 


16 3 x 48 


(2.8 fm) 3 


400 


0.176 


0.03 / 0.05 


6.81 


20 3 x 64 


(2.4 fm) 3 


559 


0.120 


0.02 / 0.05 


6.79 


20 3 x 64 


(2.4 fm) 3 


460 


0.120 


0.01 / 0.05 


6.76 


20 3 x 64 


(2.4 fm) 3 


593 


0.121 


0.007 / 0.05 


6.76 


20 3 x 64 


(2.4 fm) 3 


403 


0.121 


0.005 / 0.05 


6.76 


24 3 x 64 


(2.9 fm) 3 


136 


0.120 


0.0124/0.031 


7.11 


28 3 x 96 


(2.4 fm) 3 


261 


0.0863 


0.0062/0.031 


7.09 


28 3 x 96 


(2.4 fm) 3 


472 


0.0861 



Table 1: Ensembles used in this calculation. 

We use the Asqtad improved staggered sea quark action that has errors of 0(a s a 2 ). The im- 
proved gluon action has errors of 0{a 2 a 2 ). For the heavy valence quarks, we use the Sheikholeslami- 



Wohlert action [ ]1Q| ] (which has 0(a s a) errors) with the Fermilab interpretation To compute 
heavy quark propagators, we use point and smeared sources and sinks. The smearing approximates 
IS or 2S wavefunctions. At the sink, spatial momentum 2n / (La)[p x , p y , p z ] is given to the onium 
state. We restrict the range of p such that < 9. 

To find the onium masses, we fit two channels simultaneously for the zero momentum states. 
A delta function and a IS smearing wave function are used as the source and sink. The ground state 
and up to three excited states are included in the fit. The minimum and maximum distance from 
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the source are varied, and the best fit is selected based on the confidence level and size of error in 
the ground state and first excited state masses. After choosing the fit range, 250 bootstrap samples 
are generated to provide an error estimate. 

We must tune the hopping parameter K to the charm or bottom mass. For each lattice spacing, 
we select a sea quark mass independent value for ff. The tuning is done on an ensemble with small 
sea quark mass. In fact, as this project was done in conjuction with a study of heavy-light mesons, 
the tuning was done for the D s mass. The precision of that tuning was only about 8%. Because of 
lattice artifacts that arise for heavy states, we distinguish between the rest mass aM\ and the kinetic 
mass aMi- We use K = 0.120, 0.119 and 0.127 on the extra coarse, coarse and fine ensembles, 
respectively. The imprecision of our tuning is immediately seen in Fig. [lj 
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Figure 1: The kinetic masses of J / y/ and r\ c on each ensemble plotted as a function m q /m s the light sea to 
strange quark mass ratio. Masses are in units of r\. The physical masses are shown as lines. 



The kinetic masses have two disadvantages: their statistical errors are large compared to those 



of the rest masses, and the pattern of systematic errors is more subtle [llfl. However, for level 
splittings, a large discretization effect in the quark's rest mass drops out of the energy differences 
of hadron rest masses [|12|]. So, having tuned to approximately the right charm mass, we will now 
consider splittings based upon the rest masses of the various states. These states have been studied: 
T] c (15), 7] c (25), Y(1S), Y(2S), h c (lP), &o( lp ) and The Xci(lP) is also under study with 

a nonrelativistic P-wave source JT3[]. Currently, results for Xai^P) ai" e on ly available on one extra 
coarse ensemble. We also use the spin-averaged mass, e.g., IS = pM^n^ +M^jig\]/4 to display 
some of the splittings in the spectrum. 



3. Results 



For each lattice spacing, we plot the splittings as a function of the mass of the light sea quarks. 
A linear chiral fit is made and the splitting is extrapolated to the physical value of m = (m u + m^) /2 
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where the lattice spacing dependent value of m is determined from analysis of % and K meson 



decays constants Q14Q. The light meson decay constant analysis has not yet been completed on the 
extra coarse ensembles, so the value of m used there is only a rough estimate. 
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Figure 2: (left) The chiral extrapolation of the spin-averaged splitting between the 2S and IS states on the 
extra coarse and coarse ensembles. The extrapolated values are shown in red, and the physical value in black. 

Figure 3: (right) The splitting between the %c\ ( 3 ^i) and spin-averaged IS states on the extra coarse and 
coarse ensembles. 

Within our current statistical uncertainties, we see reasonable agreement with the experimental 
value of the splittings of the spin-averaged 2S and IS levels. The coarse value is about 2 a high. 
(See Fig. |[) As we do not yet have a full set of results for the % c2 , we cannot construct the spin 
average of the IP states. Instead, we use the % c \ and the h c . In nonrelativistic potential models, 
these two states are degenerate with each other and the spin average. The experimental splittings 
are well reproduced for these states. (See Figs. [| and f|.) 

As seen in Fig. ||J, the spin splittings are too small. For J/y and r\ c it amounts to about 10-22 
MeV. The splitting is 19% too small for the extra coarse ensemble, 14% too small for the coarse, 
and 9% too small for the fine. The splitting seems to systematically improve as the lattice spacing 
decreases. We have not yet attempted a continuum extrapolation. 

The overall agreement between this calculation with dynamical quarks and the observed spec- 
trum is very good. The most obvious issue is the smallness of spin splittings, as seen in the J/\y 
- T7 C splitting, and the mass of the %cQ state. There is some evidence of improvement as the lattice 
spacing is reduced. 



4. Outlook 

There are several ways to improve this calculation in the near future: We can include another 
fine ensemble with m q = 0.lm s . This more chiral ensemble is still being generated, but is far 
enough along that it would be worth starting the analysis. We also need to examine alternative fits 
to the ones that were selected by our automated procedure. Fermilab/MILC are almost finished 
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Figure 4: (left) Same as Fig. 3, except for the h c ( Pi). 

Figure 5: (right) Splitting between J/y/ and rj c for all three lattice spacings. There are only two ensembles 
for the fine lattice spacing, shown in purple. 
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Figure 6: Charmonium spectrum for all three lattice spacings compared with experimental values. Energy is 
offset so that zero represents the spin-averaged IS energy. From left to right for each state, crosses, octagons 
and diamonds are from the extra coarse, coarse and fine ensembles, respectively. The extra coarse %ci value 
without chiral extrapolation is the fancy cross. 



generating a new set of ensembles at a lattice spacing between extra coarse and coarse. Production 
running on additional ensembles for the new P-wave code will be done. We also plan to use 
heavier quarks to study bottomonium, which has already been studied on these configurations using 
NRQCD [g]. 

In the longer term, MILC is generating new ensembles with a m 0.06 fm that should help us 
better understand the continuum limit. However, in the current calculation, lattice spacing depen- 
dence does not seem very large compared with statistical errors for most states. 
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